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X-ray temperature measurements of clusters of galaxies are now reaching to 
redshifts high enough to constrain Q . A redshift-dependent relation that maps 
^ \ these X-ray temperatures to the virial masses of clusters is an essential ingredient 

when one is trying to determine cosmological parameters from cluster evolution. 
Most such relations assume that clusters form from top-hat perturbations that 
virialized just before the time we are observing them. The smaller Q is, the less 
accurate these relations become. Here we derive a relation between virial mass 
and cluster temperature that allows for the fact that clusters form gradually 
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and cease forming when the density of the universe drops well below the critical 
value. We show how sensitively the expected redshift distribution of clusters 
depends on the mass-temperature relation used and argue that one needs to use 



a relation that yields no evolution in the low-fio limit. 



Subject headings: cosmology: observations — galaxies: clusters: general 
X-rays: galaxies 



1. Introduction 

Studies of temperature evolution in massive clusters of galaxies are now realizing their 
promise as cosmological indicators. The most massive clusters in the universe mark the 
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rarest, highest amplitude peaks in the intial perturbation spectrum that seeded structure 
formation and provide powerful leverage on parameters such as f2 , the current matter 
density of the universe in units of the critical density p cr . If f2 = 1, the development of 
structure should have proceeded unabated from early in time to the present day. However, 
the relative lack of evolution observed in the X-ray temperature function of clusters 
indicates that the formation of structure is stalling, and thus that Q ~ 0.2 — 0.5 (Henry 
1997; Bahcall et al. 1997; Carlberg et al. 1997; Donahue et al. 1998) 

Using X-ray observations of clusters to constrain cosmological models requires 
knowledge of how the X-ray temperature of a cluster (Tx) relates to its mass. Temperatures 
of distant clusters are now straightforward to measure (e.g., Donahue 1997; Donahue et al. 
1998; Henry 1997), but relating these temperatures to cluster masses is a subtler art. In 
order to define a cluster's mass, one must decide upon an appropriate boundary separating 
the cluster itself from peripheral material not yet part of the cluster. Theorists tend to 
choose a boundary that encompasses matter that has virialized and excludes matter that is 
still falling into the cluster. The mass within this boundary is then M vir , the virial mass of 
the cluster. 

Defining a cluster's mass to be M vir has an established theoretical pedigree, but this 
definition is observationally problematic. Virial masses are a staple of the Press- Schechter 
approach to the problem of hierarchical structure formation, which has become enormously 
popular owing to its basic agreement with much more sophisticated n-body techniques (e.g., 
Lacey & Cole 1994). Press & Schechter (1974) presumed that the collapse of a spherical 
top-hat perturbation would adequately represent the formation of more complicated objects. 
In such a spherically symmetric collapse, the density of a newly formed object just after 
collapse is 187r 2 p cr ps 178p cr in an Q = 1 universe, but this contrast factor drops to as low 
as 87r 2 p cr ~ 79 p cr in an open universe with f2 <C 1 (e.g., Lacey & Cole 1993). Even if one 
could measure such density contrasts accurately, a priori knowledge of Qo would still be 
needed to determine M vir . Thus, the bridge between Tx and M vir must be built primarily 
from the theoretical side. 

Numerical models of cluster formation show that cluster temperatures should indeed be 
closely related to their masses (e.g., Evrard et al. 1996). To lowest order, cluster potentials 
that arise in such models are isothermal, with Tx oc M(r)/r (e.g., Cole & Lacey 1996). If 
the mean density within radius r is A times the critical density, then the mass within that 
radius is M(A) oc T x ^ 2 p CT 1 / 2 A -1 / 2 . Proceeding in the vein of the Press-Schechter approach, 
one can define A vir to be the density contrast of a spherical top-hat perturbation just after 
collapse and virialization, in which case M vir = M(A vir ) and T x oc M% 3 pV 3 A*/ 3 . The 
density contrast parameter A vir depends in general on fl , the redshift of virialization, and 
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the value of the cosmological constant (e.g., Eke et al. 1996; Oukbir & Blanchard 1997). 
Adopting the recent-formation approximation, which assumes that the clusters we now 
observe virialized just before their light departed, then leads to a relation like the following: 



where Q(z) is the matter density of the universe in units of p cr at redshift z, (3 is the usual 
ratio of potential depth to X-ray temperature, and h is the Hubble constant in units of 
lOOkms -1 Mpc -1 (e.g., Eke et al. 1996). In the limit of Q(z) ~ 1, we arrive at the scaling 
relation T x oc M^ 3 (l + z). 

The shortcomings of this approach are well known (e.g. Viana & Liddle 1996; 
Kitayama & Suto 1996; Eke et al. 1996). Real clusters do not form at a single moment 
but rather accrete their matter over a long period of time. If the accretion rate remains 
sufficiently high, as would be the case if Qq = 1, then the clusters we now see attained 
their observed masses relatively recently, validating the recent-formation approximation 
that led to equation (|I]). If instead fio^l, then the recent formation approximation grows 
increasingly less valid in time. Note that in the limit ri — ^ 0, cluster evolution ceases by 
z = and the M vir —T x relation should no longer depend on redshift, whereas equation (J]]) 
dictates Tx oc M^ 3 (l + z) 2 / 3 . One should also keep in mind that the temperature of a real 
cluster depends on formation epoch only insofar as the formation epoch affects the shape 
of the halo density profile and furthermore that real clusters are not necessarily isothermal. 
Comparisons of simulated clusters with equation ([J) indicate that the effective value of f3 
declines modestly but systematically with time in an open universe with Qq = 0.3, slightly 
diminishing the expected amount of temperature evolution (Eke, Navarro, & Frenk 1998). 

While the relation in equation (fj) is adequate for many purposes, one would prefer to 
use a relation valid in the low-fio limit when constraining O through cluster temperature 
evolution. One way to patch the recent-formation approximation is to define an explicit 
formation redshift Zf at which a cluster virializes. One can then integrate over the 
appropriate distribution of formation redshifts to determine the properties of observed 
clusters at redshift z (Kitayama & Suto 1996; Viana & Liddle 1996). 

Motivated by the fact that clusters form gradually, not instantaneously, we present here 
an alternative approach to characterizing the redshift evolution of the M vir — Tx relation 
which has the advantage of being qualitatively correct in the low-f2o limit. In § 2 we idealize 
clusters as forming from spherically symmetric perturbations with negative radial density 
gradients rather than from top-hat perturbations. This approximation alleviates the need 
to specify an artificial formation redshift because it explicitly accounts for how clusters 
grow with time. It also yields an analytical formula for the change in virial energy with 
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virial mass. We then use the merging-halo formalism of Lacey & Cole (1993) to derive an 
expression for the time-averaged mass accretion rate of a massive cluster, which allows us to 
determine how the M vir —Tx relation should change with time. Section 3 briefly illustrates 
some implications of this relation. 



Here we derive in three steps how the M Y - n —Tx relation should evolve. First we analyze 
how a radially stratified, spherically symmetric perturbation would virialize. Then we 
determine how the mass accretion rate of a massive cluster changes with time in a Gaussian 
density field with a power-law perturbation spectrum. Combining the results of these 
analyses, we determine how the M N1X —Tx relation evolves with time, and we normalize 
this relation using equation ([!]). For brevity and clarity, we restrict our derivation to the 
analytically tractable case of an open universe without a cosmo logical constant. 



Let us idealize the perturbation that will form a cluster as having a spherically 
symmetric density profile that declines with radius r. We then can define the Lagrangian 
coordinate M to be the mass inside a given spherical shell centered on the origin. 
The evolution of each shell depends only upon the mass inside of it. In the absence 
of a cosmological constant, these shells will obey the familar parametric solution 
r(M) = [GM(t M /27r) 2 ] 1/3 (l - cos0 M ), t = {t M /2^){e M - sin0 M ), where 9 M parametrizes 
the evolution of the shell containing mass M, which formally recollapses to the origin at 
time tM (e.g., Peebles 1993). 

Drawing an analogy between top-hat collapse and the collapse of concentric shells, we 
can approximate the virial mass of the forming cluster to be the M of the shell that has just 
recollapsed to the origin. The energy this infalling shell contributes to the total is equal to 
its potential energy at maximum expansion. Because the maximum radius of the shell that 
converges upon the origin at time t is 2[GM(t/2vr) 2 ] 1/3 , the virial energy —Eoi the cluster 
accumulates according to 



An expression for M vir (t) then provides the information needed to determine E(t) and 



2. Mass- Temperature Relation in an Open Universe 



2.1. Spherical Collapse with a Density Gradient 




(2) 



T(t) cx E(t)/M viT {t). 
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2.2. 



Mass Accretion Rate for Massive Clusters 



The merging-halo formalism of Lacey & Cole (1993) can be used to derive a mean 
mass accretion rate for massive clusters. This approach extends the Press-Schechter picture 
by considering how clusters grow via accretion of smaller virialized objects. Mass accretion 
becomes quasi-continuous for very massive clusters, and in this limit the Lacey & Cole 
formalism yields a mass accretion rate that enables us to integrate equation 

In order to analyze halo merger rates, one defines <5(x, £; M) to be the local fractional 
overdensity of the universe, smoothed on mass scale M and centered on comoving point 
x at time t. Early in time the perturbations described by S are assumed to be Gaussian, 
and in the linear regime they grow in proportion to the function D(t), which depends 
on Qq. While these perturbations remain linear, their rms amplitude on scale M can be 
expressed as a(M)D(t)/D(to). Ultimately some of them grow non-linear, and they are 
assumed to viralize when their amplitudes, extrapolated from the linear regime according 
to D(t), exceed some critical threshold S c (t). One can then trace the merger history of a 
mass parcel beginning at x from time t\ to the present by keeping track of the largest M 
for which 5(x,ti; M)D(t)/D(t\) > S c (t). This largest mass is the mass of the virialized halo 
containing the mass parcel at time t, and it jumps by an amount AM each time the halo 
merges with another of mass AM. 

Equivalently, one can analyze the growth of structure by keeping the amplitudes of 
the perturbations fixed, with rms amplitude cx(M), and tracking when they exceed the 
steadily falling threshold u(t) = 5 c (t)D(t )/D(t). In an open universe with no cosmological 
constant, Lacey & Cole (1993) show that 



where tn = (ir{l / H )(l — fi ) -3 / 2 . The parameter to thus serves as a useful surrogate for 
the time coordinate in an open universe because it asymptotically approaches a fixed value 
as growth of structure slows to a halt. 

Now we are ready to determine an expression for the halo accretion rate. Lacey & Cole 
(1993) show that the probability that the mass of a virialized halo will jump from Mi into 
an interval dM 2 at M 2 within a time interval corresponding to duj is 



W (t) = ?D(to)[l + (V*) 2/3 ] 



(3) 
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where Si = a 2 (Mi). To find the mean growth rate of a halo of mass M, we can integrate 
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over this probability distribution as follows: 

\s l ~S 2 )-^-dS l . (5) 

The variance S(M) generally declines with M, so for a rare high-mass object, S^iM^) <C uj 2 . 
In this limit, the integrand's exponential behavior restricts its main contribution to a 
narrow range where Si ~ S%, reflecting the fact that most of the accreted mass comes from 
objects much smaller than the cluster. Thus, 

c 

(6) 

If the fluctuation amplitudes obey er(M) oc M _(n+3 )/ 6 , where n is the usual power-law 
perturbation index in wavenumber space, we finally arrive at the mass accretion law 
M vir oc cu~ 3 /( n+3 ). This result is equivalent to Lacey & Cole's (1993) finding that 
M vir oc t 2 /(™+ 3 ) for tt = 1. 




2.3. Evolution and Normalization of M vir — T 



x 



The typical temperature of a cluster should be proportional to the virial energy divided 
by Myir- If we define x = [1 + {t^i/t) 2 ^] and m = 5/(n + 3), we can write M v ; r oc x _3m /5 and 



dE oc (x — l)dx m . Integrating over x then gives E oc M, 
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T x oc 
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1 + 1/m) . Thus, we obtain 



(7) 



This relation behaves properly in all the appropriate limits. The time- dependent factor 
varies like (1 + z) when Vl(z) ~ 1, like pj/ 3 oc t -2 / 3 for intermediate £l(z), and goes to a 
constant as Q(z) — > 0. To normalize this expression, we can simply match it to equation (p]) 
at early times: 



kT x = (2.76 keV)/T 
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Reducing fl at fixed M vir and t pushes T x upward because the effective formation epoch 
moves to earlier times, when the universe was denser. Note also that the implied asymptotic 
behavior of temperature evolution depends on the power-law index n. Larger amounts of 
power on small scales (larger n) lead to a less sensitive dependence of M vir on u and an 
earlier cutoff for the accumulation of mass. 
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The true test of a relation between M vir and T x is how accurately it reproduces the 
evolution of the cluster temperature function in a large numerical simulation. Showing 
that M(A) correlates well with Tx for some chosen A is not sufficient because such tests 
do not probe how A V j r changes with time. As long as clusters are nearly isothermal above 
the density contrast A, one will find Tx oc [M(A)] 2 / 3 p^/ 3 A 1 / 3 . If the goal is to constrain 
cosmological parameters by applying Press-Schechter analyses to observed clusters, then the 
M viT —T x relation one uses should yield cluster temperature functions that approximate the 
results of numerical simulations as closely as possible. To accurately reproduce behavior at 
low fl, one needs a relation like equation (|7|) that asymptotically yields no evolution. In the 
most general cases where a single value of n does not adequately represent the perturbation 
spectrum, equation (0) should still be useful as a fitting formula in which the asymptotic 
constant is a free parameter. 



3. Press-Schechter Predictions and M vir —Tx Evolution 

The generally good agreement between analytical Press-Schechter models of cluster 
evolution and numerical simulations has inspired numerous applications of the Press- 
Schechter approach to X-ray cluster surveys (e.g., Oukbir & Blanchard 1997; Henry 1997; 
Carlberg et al. 1997; Donahue et al. 1998). An accurate M wlT —Tx relation is a crucial 
element in such studies because it maps the virial masses used in the Press-Schechter 
formalism to either observed X-ray temperatures or observed X-ray luminosities (Lx) 
through an L x —T x relation. Here we will show that the consequences of using an inaccurate 
M vir —Tx relation can be quite significant. 

For the purposes of studying cluster evolution, the change in M vir —Tx with z is more 
important than its normalization. Figure 1 shows how the quantity Tx(M vir , z) /Tx{M v - iri 0) 
changes for four different relations: the one derived here for n = —1 and —2, the relation 
from equation ([!]), and the relation T x oc M^ 3 (l + z) that one obtains by assuming that the 
mean density within the virial radius is some constant multiple of the background density. 
We have also computed the mean temperatures for clusters of a fixed mass at a variety of 
redshifts by integrating over formation redshifts following Viana & Liddle's (1996) method 
and with Lacey & Cole's (1993) definition of formation epoch. The points on Figure 1 show 
how the mean temperatures of these clusters that are originally 3a perturbations evolve for 
n = — 1. Note the close agreement between these points and the n = — 1 M v - n — Tx relation 
from the present paper. 

The differences between the M vir — T x relations in Figure 1 are relatively modest in 
percentage terms, but they become greatly amplified when filtered through the exponential 
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distribution of cluster masses that emerges from Press-Schechter calculations. According to 
such calculations, the differential comoving number density of clusters dn in temperature 
interval dTx is given by 

dn /2\ 1/2 3 fioPcr(O) dlna 2 

dIV = U) 2 T x M A (T x ,z)dlnAu/ c{T: " Z) ^'"^^ ■ < 9 > 

For a power-law perturbation index n and a fiducial temperature T , one can write 
v c (Tx,z) = v c0 [u(z) /u(0)}{M viT (Tx, z)/M viT (To)]( n+3 ^ 6 , where z/ c0 or some analogous 
parameter is adjusted to fit low-redshift cluster surveys: fixing z/ c0 is equivalent to fixing 
the amplitude of the power spectrum for a given n. The derived number density of clusters 
thus depends exponentially on M vir (Tx, z). 

An M v i r — Tx relation that overestimates temperature evolution will underestimate 
the virial mass corresponding to a high-z cluster of a given temperature or luminosity. 
Because the number density of clusters above a given mass decreases with virial mass, an 
underestimate of virial mass leads one to overestimate the numbers of high- z clusters that 
ought to be in a particular temperature-limited or flux-limited sample. This effect is more 
severe in a flux-limited sample because the expected luminosity is much more sensitive to 
changes in the virial mass. 

To illustrate these effects, we have computed the surface density of clusters on the 
sky given by equation @ for Q = 0.2, n = —1, various M vir — Tx relations, the Lx—Tx 
relation from Edge & Stewart (1991), which we assume remains constant with redshift 
(e.g., Mushotzky & Scharf 1997), and a particular set of selection criteria. We restrict the 
clusters to have Tx > 5keV and Fx > 10 13 ergcm -3 s _1 within a 2.4' x 2.4' detection cell 
like that of the Einstein Extended Medium Sensitivity Survey (EMSS; Henry et al. 1992). 
We also include the effects of Galactic soft X-ray absorption, assuming a mean hydrogen 
column density 3 x 10 20 cm~ 2 . Figure 2 shows the resulting surface densities, binned in 
redshift intervals of 0.1. Note that that at z ~ 0.5 the expected cluster counts differ by a 
factor ~ 3, with more modest temperature evolution yielding fewer expected clusters. This 
factor is similar to the discrepancy Oukbir & Blanchard (1997) find between the actual 
redshift distribution of clusters in the EMSS and their expectations, suggesting that the 
entire discrepancy might stem from the M Nlx —Tx relation they use. Above z ~ 0.5, the flux 
limit becomes more restrictive than the temperature limit, and the discrepancies between 
the expected cluster counts grow more pronounced. 

Clearly, Press-Schechter approaches to modelling X-ray surveys require an accurate 
M Y - 1T —Tx relation. The relation we give here is likely to describe temperature evolution 
more accurately at late times than those based on the recent-formation approximation 
because it correctly yields no evolution as Q(z) — > 0. We are currently using this M vir —T x 
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relation to study cluster evolution in the EMSS and will report our results in a future paper. 

We would like to thank Caleb Scharf for his helpful suggestions while we were 
developing these ideas. We acknowledge NASA grants grant NAG5-3257 and NAG5-3208 
for partial support. 
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Fig. 1. — Temperature evolution for various M v - n —Tx relations. The naive relation in which 
Tx oc M^/ 3 (l + z) (long-dashed line) predicts the most evolution in Tx for a given M vir . 
Temperature evolution is more modest under the late-formation approximation in an open 
universe (dotted line) and is milder still for the relation derived in equation (|7]), given 
n — — 1 (solid line) and n = —2 (short-dashed line). The points indicate how mean cluster 
temperatures evolve in models that integrate over a range of formation redshifts Zf to find 
the temperatures of clusters at a given redshift z, given an n = — 1 perturbation spectrum. 
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Fig. 2. — Redshift distribution of clusters as predicted by different M viT —T x relations. Here 
we show expectations for clusters hotter than 5 keV with fluxes exceeding 10 13 erg cur 2 s _1 
in the EMSS detect cell. The naive relation T x oc M% 3 (1 + z) (dashed line) overpredicts 
the high- 2 clusters because it overestimates temperature evolution and maps the flux and 
temperature limits to lower values of M vir than it should. 



